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Abstract 

In this note, we derive Ewald sums for Yukawa potential for three dimensional systems 
with two dimensional periodicity. This sums are derived from the Ewald sums for Yukawa 
potentials with three dimensional periodicity [G. Salin and J.-M. Caillol, J. Chem. Phys. 
113, 10459 (2000)] by using the method proposed by Parry for the Coulomb interactions 
[D.E. Parry, Surf. Sci. 49, 433 (1975); 54, 195 (1976)]. 
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The Yukawa interaction energy between two particles is given by 

771 / \ ViVj exp(-Kr) 



(1) 



e r 



where e is the dielectric constant, k the inverse of the screening length and yi the 
"Yukawa charges" defined by the properties and the state of the system ; for instance, 
at the Debye-Hiickel approximation for electrolytes or in the Derjaguin-Landau-Verwey- 
Overbeek (DLVO) theory of colloids, k and yi are related to physical parameters of systems 
as 



where a, is the diameter of the hard core of the ions, in the Debye-Hiickel approximations, 
or the radius of macroions, in DLVO theory, and p and q are respectively the density of 
ions or counterions and their charge, ks the Boltzmann constant and T the temperature. 
Yukawa interactions between particles are used in numerical simulations as effective po- 
tentials to simulate systems as plasmas, dusty plasma, colloids, etc. ; on general ground, 
such potentials may be used as a reasonable approximation, as soon as some micro- 
scopic degrees of freedom may be approximated to a continuous background leading to a 
screening of the direct interaction between particles, while the spherical symmetry of the 
interaction is preserved. 

As outlined in ref . [1] , if k is large enough, the screening length can be much smaller than 
simulation box lengths, then interactions between particles are not long ranged and, in 
practice, a simple truncation of the potential, with the use of the minimum image conven- 
tion, could be sufficient. On the contrary, if k is not large or quite small, then interactions 
between particles may be long ranged and images of particles introduced by the periodic 
boundary conditions may contribute significantly to the energy of the system. In these 
cases, a crude truncation of the potential could lead to strong bias in computations (for 
Coulomb interactions, see for instance refs. [2-4] for errors introduced by crude truncations 
of long ranged potentials). To handle these latter cases, an Ewald method for systems with 
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three dimensional periodicity and Yukawa interaction potentials has been exhibited 1 . 
Many interesting systems which interaction between particles can be approximated by 
Yukawa potentials are also confined to quasi-two dimensional geometries 5 ' 6,7 , therefore 
an Ewald method is of interest to permit to simulate the properties of these quasi-two 
dimensional systems for any value of the k parameter including at low counterions con- 
centration or high temperatures. 

In this note, we derive Ewald sums for Yukawa potential in quasi-two dimensional sys- 
tems from results of ref.[l] following the same derivation done by Parry 8 for Coulomb 
interactions. For Coulomb interaction in quasi-two dimensional systems several methods 
exist 9 ' 10 , in particular some methods used the Ewald method for three dimensional sys- 
tems with a highly asymmetric box 11 ' 12 and by adding correction terms related to the 
total dipole of the simulation box ; a general review on Coulomb interaction in quasi 
two-dimensional systems is done in ref.[10]. In a forthcoming work, some numerical im- 
plementations on a test system will be given ; the present work is devoted only to provide 
a simple derivation of Ewald sums for Yukawa potential in quasi-two dimensional systems. 
As computed by Salin and Caillol 1 , the Ewald- Yukawa interaction energy is given by 

E = E r + E k - E Se[i (2) 

with the short ranged contribution 

1 ^ ' D(r ij} n,K;a) 

Er = -Y.Y.y^ lr .. + nLl ( 3 ) 

where 

D(rij, n, k\ a) = erfc(a | + nL | H ) exp(/t | + nL |) 

2a ( 4 ) 
+erfc(o; | + nL \ — —— ) exp(— k, \ + nL |) 

and the distance between the pair (i, j) of particles, the long ranged contribution 

2tt _ exp(-fc 2 + /t 2 )/4a 2 ),_ |2 
E k = y22 -£T— 2 |2>exp(*fc.r,)| (5) 
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and the self interaction 



E Seli 



a 



-p(-^)-2 erfc (^)jE^ 



(6) 



where we have set e = 1 and used conventional notations for Ewald sums ; namely, V is the 
volume of the simulation box, nL the condensed notation for the vectors of the periodic 
boundary conditions, k are the vectors belonging to the reciprocal lattice associated with 
the three dimensional periodicity and a the damping parameter of the Ewald method. In 
Eq.(3), the prime in the sum over n indicates that for n = 0, the self terms i = j are not 
included. The short ranged contributions for the Ewald sums for Yukawa potentials in 
periodic three dimensional systems are given as in Eq.(3), for quasi-two dimensional 8 the 
short ranged contributions are also given by Eq.(3) ; for pratical applications, a choice of 
the Ewald damping parameter a is done such that summations over images are restricted 
to the minimum image convention 1 . 

In the work by Parry 8 , the Ewald sums for quasi-two dimensional systems are derived 
from the Ewald method for three dimensional systems by letting the spatial periodicity 
along the third direction to tend to infinity (L z — > oo). 

In the following, we use the notations k = G + ke z where G are the vectors belonging to 
the reciprocal lattice associated with the two dimensional periodicity and k = 2nm/L z , 
(m integer), where L z is the spatial periodicity of the simulation box along the e z ; we set 
also Vij = Sij + Zije z and V = AL Z , where Sij is the distance in the plane perpendicular 
to e z and A is the surface of the simulation box for quasi-two dimensional systems. 
Following Parry 8 , the long ranged contribution is separated into two contributions. The 
first contribution 81 ^, noted Eq^ , is obtained for G^O and the second 8 '- 6 -' for G = 0, as 
a reminiscent contribution of summation over k as L z — > oo, noted below -^^=0- 
From Eq.(5), we have 




YjViVj Yl exp(iG.s ii )exp(-(G 2 + n 2 )/Aa 2 ) 
v G^o 
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Taking into account the identity 13 
1 +2? exp(-A; 2 /4a^ 



we found 



1 f + °° exp(ikz) exp(— k 2 /4a 2 ) 

d fx ~ 



[ 

J — < 



k 2 + (G + k 2 



E, 



GMo *J 



(8) 



(9) 



with 



1 / \ 
F(K,z;a) = — [exp(Kz)erfc( \- az) + exp(—Kz)eric( az)j 

The second contribution is given by 



(10) 



E, 



(b) 
G=0 



— Y,yiyjF(K, Zij ;a) 



(11) 



As it is shown explicitly on Eqs.(9) and (11), the long ranged contributions of Ewald- 
Yukawa sums for quasi-two dimensional systems can not be reduced to one particle sum- 
mation, because of the complicated dependence on z^. 
Self energy is given by 



E Q2D {Se\f) = 



-= exp(— K 2 /4a 2 ) — — erfc(/t/2o;) 

V 2 



E^ 2 



(12) 



7r erfc(/t/2o;) ^ , 
£ G=0 (Self) =- y -^Y,yt 



A K 

L 

The Ewald sums for Yukawa potential in quasi-two dimensional systems are given by 



E = E r + 4 a io + 4=o - £ Q2D (Self) - E G=Q (Se\i) 



(13) 



with each contribution given by Eqs. (3,9,11) and (12). 

When k — > 0, the Helmholtz equation, which solutions lead to Yukawa potentials, becomes 
a Poisson equation, which solutions lead to Coulomb potentials 1 ; then the Ewald sums for 
Coulomb interactions may be recovered from Eqs. (3,9,11) and (12), under the condition 
that the electroneutrality of the periodic system is recovered too 14 ' 15 (for instance by 



letting the counterions concentration acting as an uniform continuous background). 
When Zij = for all pairs of particles, then all particles are confined in a plan. In this 
limit, the Ewald sums for a quasi-two dimensional systems reduce to the Ewald sums of 
two dimensional systems. Then, Eqs.(9,12) and (13) become 



hi — his + hiQ^Q — hi, 



Self 



(14) 



with 



E, 



^Self 



A G*o VG 2 + k 2 



a 2 I a 2\ K 



—^=exp(—K /4a ) — — erfc(«;/2a) 



(15) 



and E s given by Eq.(3), with Sij instead of r^. 

These results may be applied directly in computer simulations, nevertheless the efficiency 
suffers from the lack of reduction in one particle summations, as for Coulomb potentials 
in quasi-two dimensional systems, because of the complicated dependence on z^. For 
molecular dynamics implementations, the forces acting on particles can easily be obtained 
with computations of the derivatives of the interaction potential. 
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